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Mathematical mean-field approaches play an important role in 
different fields of Physics and Chemistry, but have found in recent 
works also their application in Economics, Finance and Game The- 
Q^ ' ory. The objective of our paper is to investigate a special mean-field 

Q, , problem in a purely stochastic approach: for the solution (Y, Z) of 

a mean-field backward stochastic differential equation driven by a 
forward stochastic differential of McKean-Vlasov type with solution 
Cu ■ Jf we study a special approximation by the solution (X'^ ,Y'^ ,Z'^) 

of some decoupled forward-backward equation which coefficients are 
governed by A'' independent copies of {X^ ,Y^ ,Z^). We show that 
the convergence speed of this approximation is of order 1/VN. More- 
over, our special choice of the approximation allows to characterize 
the limit behavior of ^/]V(X'^ -X,Y'^ - Y, Z'^ - Z). We prove that 
this triplet converges in law to the solution of some forward-backward 
stochastic differential equation of mean-field type, which is not only 
f^ ■ governed by a Brownian motion but also by an independent Gaussian 

field. 
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t^^ . 1. Introduction. Our present work on a stochastic limit approach to a 
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mean-field problem is motivated on the one hand by classical mean-field ap- 
proaches in Statistical Mechanics and Physics (e.g., the derivation of Boltz- 
mann or Vlasov equations in the kinetic gas theory), by similar approaches 



?H ' in Quantum Mechanics and Quantum Chemistry (in particular the density 
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functional models or Hartree and Hartree-Fock type models) but also by 
a recent series of papers by Lasry and Lions (see [7] and the references 
therein). On the other hand it stems its motivation from partial differen- 
tial equations of McKean-Vlasov type, which have found a great interest in 
recent years and have been studied with the help of stochastic methods by 
several authors. 

Lasry and Lions introduced in recent papers a general mathematical mod- 
eling approach for high-dimensional systems of evolution equations corre- 
sponding to a large number of "agents" (or "particles"). They extended the 
field of such mean-field approaches also to problems in Economics and Fi- 
nance as well as to the game theory: here they studied A'^-players stochastic 
differential games [7] and the related problem of the existence of Nash equi- 
librium points, and by letting N tend to infinity they derived in a periodic 
setting the mean-field limit equation. 

On the other hand, in the last years models of large stochastic particle 
systems with mean-field interaction have been studied by many authors; 
they have described them by characterizing their asymptotic behavior when 
the size of the system becomes very large. The reader is referred, for ex- 
ample, to the works by Bossy [1], Bossy and Talay [2], Chan [4], Kotelenez 
[6], Meleard [8], Overbeck [10], Pra and Hollander [14], Sznitman [15, 16], 
Talay and Vaillant [17] and all the references therein. They have shown that 
probabilistic methods allow to study the solution of the linear McKean- 
Vlasov PDE (see, e.g., Meleard [8]). With the objective to give a stochastic 
interpretation to such semilinear PDEs we introduced in [3] the notion of 
a mean-field backward stochastic differential equation: backward stochastic 
differential equations (BSDEs) have been introduced by Pardoux and Peng 
in their pioneering papers [11] and [12], in which they proved in particular 
that, in a Markovian framework, the solution of a BSDE describes the vis- 
cosity solution of the associated semilinear PDE. The reader interested in 
more details is also referred to the paper by El Karoui, Peng and Quenez [5] 
and the references therein. However, in order to generalize this stochastic 
interpretation to semilinear McKean-Vlasov PDEs we have had to study a 
new type of BSDE, which takes into account the specific, nonlocal struc- 
ture of these PDEs; we have called this new type of backward equations 
mean-field BSDE (MFBSDE). 

The objective of the present paper is to investigate a special mean-field 
problem in a purely stochastic approach: for the solution (Y, Z) of a MF- 
BSDE driven by a forward SDE of McKean-Vlasov type with solution X 
we study a special approximation by the solution {X'^ ,Y^ ,Z^) of some 
decoupled forward-backward equation which coefficients are governed by A'^ 
independent copies of {X ,Y , Z ). We show that the convergence speed 
of this approximation is of order l/\/iV. Moreover, our special choice of 
the approximation allows to characterize the limit behavior of \/N{X — 
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X, Y — Y, Z — Z) . We prove that this triplet converges in law to the so- 
lution of some forward-backward SDE of mean-field type, which is not only 
governed by a Brownian motion but also by an independent Gaussian field. 
To be more precise, for a given finite time horizon T > 0, a d-dimensional 
Brownian motion W = (M^t)tG[o,T] and a driving d-dimensional adapted 
stochastic process X = (^i)te[o,T] we consider the BSDE of mean-field type 

dYt = -E[f{t,X,At)U=A,dt + ZtdWt, t€[0,T], 
(1.1) 

YT = EMx,XT)h=Xr 

with At = {Xt,Yt,Zt), where {Y,Z) = (^t, ^t)tg[o,T] denotes the solution of 
the above equation. This equation also can be written in the form 

dYt = - f f{t,At,\)PA,{dX)dt + ZtdWt, te[0,T], 

Yt= I ^XT,x)]PxAdx). 

Such type of BSDEs has been studied recently by Buckdahn, Li and Peng 
[3] in a more general framework. In [3] for the "Markovian-like" case in 
which the process X is the solution of a forward equation of McKean-Vlasov 
type (of course, X is not a Markov process) it has been shown that such 
a BSDE gives a stochastic interpretation to the associated nonlocal partial 
differential equations. In the present work, given an arbitrary sequence of 
adapted processes X^ = (^j^)fe[o,T] such that -E[suptg[o,T] \^t^ ~ ^tP] ~^ Oj 
we study the approximation of the above MFBSDE by a backward equation 
of the form 

N 



(1.2) 



N 



hT.'^i^' 






This approximating BSDE is driven by the i.i.d. sequence A^'-^ = {X^'^ ,Y^'^ , 
Z^'^), 1 < j < N, of triplets of stochastic processes, which are supposed to 
obey the same law as A^ = (X^,y^,Z^) and to be independent of the 
Brownian motion W. In Section 3 we will introduce a natural framework in 
which such BSDEs can be solved easily, and we will show that the solution 
{Y , Z ) converges to that of (1.1) (Theorem 3.1). For the study of the 
speed of the convergence we assume that the process X is the solution of 
the stochastic differential equation (SDE) of McKean-Vlasov type 

dXt = E[a{t,x,Xt%=XtdWt + E[b{t,x,Xt%=Xtdt, 
(1.3) 

Xo = xo£R"', t£[0,T], 



(1.5) 
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and we approximate X in the same spirit as (Y, Z) by the solution X of 
the forward equation 

dX!" = ^ E ^(*' ^f > ^f ") dWt + j^Y: bit, ^t\ X^'^) dt, 

(1.4) 

X^ = xo, te[0,T]. 

We point out the conceptional difference between our approximation and 
the standard approximation for the McKean-Vlasov SDE which considers 
as X the dynamics X '^ of the first element of a system of N particles 
which dynamics is described by the system of SDEs 

dX^' = ^ E -(^' ^r^ ^f '^■) dWi + 1 E bit, Xf ^ Xf '^■) dt, 

X^'* = xo, tG[0,T], l<i<iV, 

governed by the independent Brownian motions W^, . . . , W^ (see, e.g.. Bossy [1], 
Bossy and Talay [2], Chan [4], Meleard [8], Overbeck [10], Pra and Hollander 
[14], Sznitman [15, 16], Talay and Vaillant [17]). 

By proving that the speed of the convergence of A^ = (X^, 1"^, Z^) to 
A = {X,Y,Z) is of order 1/yN (see Proposition 3.2 and Theorem 3.2) we 
know that the sequence y/N{X^ - X,Y^ -Y,Z^ - Z), N >1, is bounded, 

sup (|Xf - Xt\^ + \Yt^ - Yt\^) + / |Zf - Ztl^dt] < +CX), 
Ue[o,T] Jo i 

so that the question of the limit behavior of the sequence \^{X^ — X, Y^ — 
Y, Z^ — Z), N > 1, arises. Our special choice of the approximating se- 
quence turns out to be particularly helpful for the study of this question. 
We show that the sequence ^/N{X^ — X, Y'^ — Y, Z^ — Z) converges in law 
on_the space C{[0,T];R'^) x C([0,r];M) x L^[0,T];R'^) to the unique solution 
(X,Y,Z) of a (decoupled) forward-backward SDE of mean-field type 

Xt= f\i'\X,)ds+ j\f\X,)dWs 

JO JO 

(1.6) + f {E[{VMx,Xs)]a^=x.Xs + E[{V^,h){x,Xs)Xs\x=x.)ds 

Jo 

+ / {E[{V,a){x,X,)l=XsXs + E[{V^>a){x,Xs)Xsh=Xs)dWs, 
Jo 

Yt = {C^^\Xt) + E[V^^X,XT)U=XrXT + E[V^,Hx,XT)XTl=Xr} 
T 



sup E 

N>1 



+ J^ {^i:'HAs)+E[Vxf{X,As)]x=AAs 
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(1-7) _ 

+ E[Vx'f{X, As) As]x=A J ds 

which is not only driven by a d-diniensional Brownian motion W but also by 
a zero-mean (2d + l)-parameter Gaussian field ^ = {Q {x),S,t {x),^^^^{x), 
Q ' {x, y, z), {t, X, y, z) e [0, T] x M"^ x M x W^} which is independent of W and 
whose covariance function is defined in Theorem 4.2. 

Such a characterization, and in particular the study of convergence in law 
of a sequence of solutions of BSDEs to the solution of another BSDE, which 
includes not only the y-component but also the Z-component and is not 
embedded in a Markovian framework is new to our best knowledge (observe 
that the solutions of SDEs of mean-field type are not Markovian). The main 
idea for the proof of the convergence result consists in rewriting the forward- 
backward equations for A^ = {X^ ,Y^ ,Z^) by using the stochastic field 

e^ = (^^'^,e^'^,e^'^,e^'^), 

el'^(x) = ^(^(x,xf '^^) - £;[7(x,xf '^)]), (t,x) G [o,r] X M-^ 

with 7 = 6 for i = 1, and 7 = 0" for i = 2, 

1 ^ 
e^^{x) = ^=Y.{^{x,X^^') - E[^{x,X^^')]), X G M^ 

1 ^ 
^t'^(A) = 4? E(/(A'^^'') - ^[/(A,Af '')]), (t, A) G [0,r] X R2<i+\ 

and in proving that this stochastic field converges in law to ^. The proof 
of this convergence in law on C([0,r] x M^'^+-^;M'^+'^ +^) (Proposition 4.3) 
is split into the proof of the tightness and that of the convergence of the 
finite-dimensional laws. For getting the tightness of the laws of i ,N > 1, 
we have to suppose that the function f{{x,y,z), {x',y',z')) does not depend 
on z' , that is, there is no averaging with respect to the Z-component of the 
BSDE; in the proof of the convergence of the finite-dimensional laws the 
central limit theorem plays a crucial role. Once having the convergence in 
law of ^ to ^ we can apply Skorohod's representation theorem in order 
to obtain the almost sure convergence ^'^ — > ^' for copies of ^^ and ^ on 
an appropriate probability space. This allows to show that the associated 
redefinitions of the triplets y/N{X^ - X,Y'^ -Y,Z'^ - Z), N > 1, converge 
on this new probability space to the solution {X , Y, Z) of the redefinition of 
the system (1.6)-(1.7) (Proposition 4.4). Our main result follows then easily. 
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Our paper is organized as follows: the short Section 2 recalls briefly some 
elements of the theory of backward SDEs which will be needed in what 
follows. In Section 3 we introduce the notion of MFBSDEs and the frame- 
work in which we study them, and we prove the existence and uniqueness. 
Moreover, the approximation of MFBSDE (1.1) by (1.2) is studied, and 
in a "Mar kovian- like" framework in which (1.1) and (1.2) are associated 
with the forward (1.3) and (1.4), respectively, the convergence speed is es- 
timated. Finally, Section 4 is devoted to the study of the limit behavior of 
the triplet y/N{X^ - X,Y'^ -Y,Z^ - Z), N >l. In order to give to our 
approach a better readability we discuss first the limit behavior of the triplet 
\/N{X — X), N >1, to investigate later the limit behavior of the triplets 
by using analogies of the arguments. 

2. Preliminaries. The purpose of this section is to introduce some basic 
notions and results concerning BSDEs, which will be needed in the subse- 
quent sections. In all that follows we will work on a slight extension of the 
classical Wiener space {Q,J^,P): 

• For an arbitrarily fixed time horizon T > and a countable index set / 
(which will be specified later), Q is the set of all families {uj^)i^j of contin- 
uous functions w* : [0,r] -^ W^ with initial value (O = Co([0,r];M°')^); it is 
endowed with the product topology generated by the uniform convergence 
on its components Co([0,T];M'^); 

• ;B(r2) denotes the Borel cr-field over 0, and B = (VF*)jg/ is the coordinate 
process over J7: W/(a;) =ujl,te [0,T],uj e^,i ^ I; 

• P is the Wiener measure over (il., B{^)), that is, the unique probability 
measure with respect to which the coordinates W^,i € /, form a family of 
independent ti-dimensional Brownian motions. Finally, 

• ^ is the (7-field B{Q) completed with respect to the Wiener measure P. 
Let W := W^. We endow our probability space {^,J-,P) with the filtration 
F = {J^t)tG[o,T] which is generated by the Brownian motion W, enlarged by 
the (7-field Q = a{W},t £ [0,T],i £ I\ {0}} and completed by the collection 
Mp of all P-nuU sets: 

Tt = J'}^vg, te[0,T], 

where F^ = {T}^ = a{Wr,r < t} VAfp)telo,T] ■ We observe that the Brownian 
motion W has with respect to the filtration F the martingale representa- 
tion property, that is, for every .T^j'-measurable, square integrable random 
variable ^ there is some d-dimensional F-progressively measurable, square 
integrable process Z = (.^f)tg[o,T] such that 

C = E[C\g]+ f ZtdWt, P-a.s. 
Jo 
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We also shall introduce the following spaces of processes which will be 
used frequently in the sequel: 



S'p([0,T]) = < (^t)tg[o,T] continuous adapted process: 



E 



sup \Yt\'^ 
.te[o,T] 



< +CX) 



Lp([0,T];]R ) = \ (■^t)te[o,T]I^ -valued progressively measurable process: 



E 



T 

\Zt\'^dt 





< +00 



(Recall that \z\ denotes the Euclidean norm of z G M".) Let us now consider 
a measurable function g:Q x [0,T] x M x M*^ — > R with the property that 
((^(t,?/, 2:))tg[o,T] is F-progressively measurable for all {y,z) in M x M.'^. We 
make the following standard assumptions on the coefficient g: 

(Al) There is some real C > such that, P-a.s., for all t G [0,T], 2/1,2/2 G 

\g{t,yi,zi) -g{t,y2,Z2)\ < C(|yi -2/2I + 1^1 -Z2I). 

(A2) g{;0,0)£LU[0,T];R). 

The following result on BSDEs is by now well known; for its proof the reader 
is referred, for instance, to the pioneering work by Pardoux and Peng [11], 
but also to El Karoui, Peng and Quenez [5]. 

Lemma 2.1. Let the coefficient g satisfy the assumptions (Al) and (A2). 
Then, for any random variable S, G L'^{^,J-t-,P), the BSDE associated with 
the data couple {g,^) 

Yt = C + J g{s,Ys,Zs)ds- J ZsdBs, 0<t<T, 

has a unique F-progressively measurable solution 

(y,Z)GV([0,T])xL|([0,r];M'^). 

The proof of Lemma 2.1 is related with the following, by now standard 
estimate for BSDEs. 

Lemma 2.2. Let (51,^1), (325^2) be two data couples for which we sup- 
pose that gk satisfies the assumptions (A.l) and (A. 2) and^k £ L'^{^)^t-,P): 
k = 1,2. We denote by {Y , Z ) the unique solution of the BSDE with the 
data {gk,ik)^k = 1)2. Then, for every 5 > Q, there exist some 7(= 75) > 
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and some C{= Cs) > only depending on 5 and on the Lipschitz constants 
of gkyf^ = ^i'^j such that, with the notation 

{Y,'Z) = {Y^-Y\Z^-Z^), 9 = 91-92, ^ = 6-6, 

we have 

e\ r e-y\\Y,\^ + \Zt\^)dt]<CE[e^^m+6E\ r e^'\g{t,Y,\zl)\'dt 
Jo J Uo 



Besides the existence and uniqueness result we shall also recall the com- 
parison theorem for BSDEs (see Theorem 2.2 in El Karoui, Peng and Quenez 
[5] or also Proposition 2.4 in Peng [13]). 

Lemma 2.3 (Comparison theorem). Given two coefficients gi and g2 
satisfying (Al) and (A2) and two terminal values 1^1,^2 G L^{^7^T,P), we 
denote by {Y^,Z^) and (y^,Z^) the solution of the BSDE with the data 
(^1,(71) and {^2-, 92), respectively. Then we have: 

(i) (Monotonicity) . If ii>^2 o-nd gi > g2,a.s., then Y^ > Y^ , for all 
te [0,T], a.s. 

(ii) (Strict monotonicity). If, in addition to (i), also P{6.i > £,2} > 0, then 
we have P{Y^^ > Y^] > {), for all Q <t < T, and in particular, Y^^ > Yf^. 

After this short and very basic recall on BSDEs let us now investigate the 
limit approach for mean- field BSDEs (MFBSDEs). 

3. Mean-field BSDEs. 

3.1. The notion of mean-field BSDEs. Existence and uniqueness. The 
objective of our paper is to discuss a special mean-field problem in a purely 
stochastic approach. Let us first introduce the framework in which we want 
to study the limit approach for MFBSDEs. For this we specify the countable 
index set introduced in the proceeding section as follows: 

/ := {i\i G {1, 2, 3, . . .}^, A; > 1} U {0}. 

For two elements i = {ii, . . . ,ik), i' = {i'l, • • • ,i'^,/) of I we define i ® i' = 
{ii, . . . ,ik,i'i, . . . ,i'j^,) £ I, [with the convention that (0) ® z = i]. Then, in 
particular, for all i>l, (i) (Bi = {i, ii, • • • , ik)- 

We also shall introduce a family of shift operators Q^:Q -^ Q, k>0, 
over 0. For this end we set Q^{uj) = {uj^^^®^)j^i,uj G $7,^ > 0, and we ob- 
serve that Q^{lo) can be regarded as an element of $7 and 0^ as an operator 
mapping O into Q. The fact that all these operators Q -.Q — > O let the 
Wiener measure P invariant (i.e., P@fe = P) allows to interpret B'^ as oper- 
ator defined over L°{n,T,P): putting e'=(0(w) := e(0''(w)),cj G ft, for the 
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random variables of the form ^(lj) = /(tJ^J, . . . ,W("), ii, . . . ,in S /,ii, • • • ,tn S 
[0,T],/ G C{R'^^"'),n > 1, we can extend this definition from this class of 
continuous Wiener functionals to the whole space L^{Q,J-,P) by using the 
density of the class of smooth Wiener functionals in LP{U.,J^,P). We ob- 
serve that, for all ^ G L^{^,J^,P), the random variables @''{^),k > 1, are 
independent and identically distributed (i.i.d.), of the same law as ^ and 
independent of the Brownian motion W. 

Finally, to shorten notation we introduce the {N + l)-dimensional shift 
operator Gat = (0*^,0^, .. .,0^), which associates a random variable ^ G 
L^{Q,J^,P) with the (A^ + l)-dimensional random vector @NiC) = (^,©"'^(0) 
...,0^(^)) (notice that 0*^ is the identical operator). If ^ on its part is a 
random vector, 0^(^) and &n{0 ^-^^ defined by a componentwise application 
of the corresponding operators. 

For an arbitrarily fixed natural number A^ > let us now consider a 
measurable function /:[7 x [0,r] x M^+^ x m(^+i)x"' -^ R with the prop- 
erty that (/(i,y,z))(g[o^T] is F-progressively measurable for all (y,z) in 
]^7V+i ^ ^(Af+i)x(i_ -^g make the following standard assumptions on the co- 
efficient /, which extend (Al) and (A2) in a natural way: 

(Bl) There is some constant C > such that, P-a.s., for all t G [0,T],yi, 
y2GM^+\zi,Z2GM(^+^)^^ 

|/(t,yi,zi)-/(t,y2,Z2)| <C(|yi-y2| + |zi -Z2I). 

(B2) /(•,0,0)GL2(0,r;M). 

The following proposition extends Lemma 2.1 to the type of backward 
equations which will be used for the approximation of MFBSDEs. 

Proposition 3.1. Let the function f satisfy the above assumptions 
(Bl) and (B2). Then, for any random variable ^ G L'^{^,J-'t,P), the BSDE 
associated with (/, ^) 

(3.1) dYt = -f{t,QN{YuZt))dt + ZtdWt, tG[0,r], Yt = i 

has a unique adapted solution {Y,Z) G 5'|([0,r];]R) x L|,([0,r];M'^). 

Proof. Let H"^ := Ll{[0,T];R) x L|([0,r];IR'^). It is sufficient to prove 
the existence and the uniqueness for the above BSDE in H^. Indeed, if 
(Y, Z) is a solution of our BSDE in H^ an easy standard argument shows 
that it is also in B^ := 5|([0,r];M) x L|([0,r];M'^). On the other hand, the 
uniqueness in H^ implies obviously that in its subspace B^. 

For proving the existence and uniqueness in H^ we consider for an arbi- 
trarily given couple of processes {U, V) G H^ the coefficient gt = f{t, Q^iUt, 
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Vt)), t G [0,T]. Since g is an element of Lp([0,T];M) it follows from Lemma 
2.1 that there is a unique solution <^{U, V) := {Y, Z) £ H^ of the BSDE 

dYt = -f{t,eN{Ut,Vt))dt + ZtdWt, tG[0,T], Yt = C. 

For a such defined mapping $ : H^ -^ H^ it suffices to prove that it is a 
contraction with respect to an appropriate equivalent norm on H^ in order to 
complete the proof. For this end we consider two couples {U^,V^), {U'^, V"^) £ 
H and (F^, Z'') = ^{U'',V''),k = 1, 2. Then, due to Lemma 2.2, for all 6 > 
there is some constant 7 > (only depending on 6) such that, with the 
notation (Y,Z) = {Y^ -Y^,Z^- Z^), 



E 



re"'\\Yt\^ + \Zt\^)dt 
Jo 



<5E 



e^yit,eMiul,v,')) - fit,eMiulVt^))\''dt 



Let (17, F) = {U^ - U^,V^ - V^). Then, from the Lipschitz continuity of 
f{uj,t,-,-) [with Lipschitz constant L which does not depend on (u!,t)] and 
the fact that the random vectors G'=(I7,F) (= e'=([/\ V^) - &''{U'^,V^)),0 < 
k < N obey all the same law, we have 

2 



E 



re^'\\Yt\^ + \Zt\^)dt\ <5L^E [ 
JO . Jo 



N 



dt 



k=0 

rp TV 

<5L\N + l)f e^'y^E[\Q\UuVt)\^]dt 
^0 k=0 

= 6L\N + lf f e^'E[\Ut\'^ + \Vt\^]dt 
Jo 

Jo 



t\^)dt 



for 6 := 2L "^{N + 1) ^. This shows that if we endow the space H^ with the 
norm 



11(^,^^)11 



m 



E 



re^\\Ut\' + \Vt\') 
Jo 



dt 



1/2 



{U,V)£H^, 



the mapping ^-.H^ — > H^ becomes a contraction. Thus, the proof is com- 
plete. D 

We now introduce the framework for the study of the limit of the above 
BSDE as A'' tends to +00. For this end let be given a data triplet ($,5, A') 
with the following properties (CI), (C2) and (C3): 
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(CI) 5 : i7 X [0, T] X (M™ x E x R'^)^ ^ R is a bounded measurable function 
which is Lipschitz in (u, v) e (M"^ x R x R*^)^ with a Lipschitz constant C, 
that is, P-a.s., for all t G [0,r] and (u, v), (u', v') G (R'" x R x R'^)^, 

|5(t,(u,v))-<7(t,(u',v'))|<C(|u-u'| + |v-v'|); 

(C2) <1> : il X R™ X R™ ^ R is a bounded measurable function such that 
^{lo,-,-) is Lipschitz with a Lipschitz constant C, that is, P-a.s., for all 

(x,f), {x',x')£R'^, 

\^{x,x) -$(x',x')| <C(|x-x'| + \x-x'\). 

(C3) X = (X^)7v>i is a Cauchy sequence in 5|([0,r];R™), that is, there 
is a (unique) process X G 5p([0,T];R'") such that 



E 



sup |Xi^-Xt|2 

t6[0,T] 



as7V^+oo. 



Remark 3.1. In the next section we will consider as sequence X^ , 
N >1, a special approximation of the solution of a forward SDE of McKean- 
Vlasov type. This special choice will allow to study the convergence speed 
of the BSDEs as N tends to +00 and to characterize the nature of this 
convergence more precisely. 

Given such a triplet {^,g,X) satisfying the assumptions (C1)-(C3) we 
put 

1 ^ 
f''{u;,t,y,z):=-Y,g{Q''{u;),t,Xf{u;),{yo,zo),X,''{eHco)),{yk,Zk)) 

fc=i 

for {io,t)enx [0,T],y = (yo,...,yiv)GM^+\z = (zo,...,^7v)eK(^+')x^ 
and 

1 ^ 

^'^H^=^E^(0'M,^tH,^t(0'M)), n>i. 
k=i 

We observe that, for every A^ > 1, /^ satisfies the assumptions (Bl) and 
(B2). Thus, due to Proposition 3.1, for all A^ > 1, there is a unique solution 
(Y^,Z^) of the BSDE(iV) 

(3.2) Yf = C^ + J^^ f^is,@N{Yl',Z^))ds- 1^^ Z^ dWs, tG[0,r]. 

We remark in particular that the driving coefficient of the above BSDE(A^) 
(3.2) can be written as follows: 

f^{s,QNiYl',Z^)) = ^J:ie''g)is,{X^,Yf,Z^),eHx^,Yj',Z^)), 

k=i 
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sG [0,r] [recall that Q''g{oj,s,{u,v)) := g{@''{uj),s,{vL,v))]. Our objective 
is to show that the unique solution of BSDE(A/') (3.2) converges in B^ = 
5|,([0,r]) X L|([0,r];]R'') to the unique solution {Y,Z) of the MFBSDE 



(3.3) 



dYt = -E[g{t,u,At)]\^=^,dt + ZtdWu te[0,T], 



where we have used the notation A = (X, Y, Z). 



Lemma 3.1. Under the assumptions (C1)-(C3) the above MFBSDE 
(3.3) possesses a unique solution {Y,Z) G B^. 

Since the proof is straightforward and uses essentially the argument de- 
veloped in the proof of Proposition 3.1, we omit it. See also [3]. 
We now can formulate the following theorem: 



Theorem 3.1. Under the assumptions (C1)-(C3) the unique solution 
{Y^,Z^) of BSDE(N) (3.2) converges in B^ to the unique solution {Y,Z) 
of the above MFBSDE (3.3): 



E 



sup \Y^ -Yt\ + / \Z^ - Zt\ dt 
;p[n_Ti JO 



Ue[o,r] "'0 

Proof. We first notice that 







as N ^ +CXO. 



Step 1. For all p>2, 



E 



E 



1 ^ 
-Y,{e''g){t,At,e'{At)) - E[g{t,u,At)]\^^ 

k=i 

1 ^ 

-Y.ie'^<^){XT,eHxT))-EMx,XT)]\.=xr 
k=i 



dt 



PI 







as N ^ +00 [recaU that (e'=$)(a;,XT,e'=(XT))(w) := $(e'=(w),XT(w), 

XrieH^m- 

For proving the first convergence we consider arbitrary t G [0,T] and u G 
R*" X M X M . Observing that the sequence of random variables {& g) {t, u, Aj) , 
k>l, is i.i.d. and of the same law as g{t,u,At), we obtain from the Strong 
Law of Large Numbers that 

N 



l-J2i&'9){t,u,At)^E[g{t,u,At)], 



N 



fc=i 
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P-a.s., as A^ — > +00. Let now, for an arbitrarily small e > 0, Af : J7 ^ M™" x 
M X M be a random vector which has only a countable number of values 
and is such that |At — A|| <e, everywhere on Q. Then, obviously, 

1 ^ 



N 



fc=i 



P-a.s., as A^ tends to +00, and from the Lipschitz continuity of g(w,i, •, v), 
which is uniform in {u!,t,v), it follows that we also have the convergence for 
At instead of Af : 

1 ^ 

k=i 

P-a.s., as A^ ^ +00. Finally, in view of the boundedness of g and, hence, 
of that of the convergence, we get the announced result. With the same 
argument we also get the L^-convergence for the terminal condition, for all 
p>2. 

Step 2. By applying Lemma 2.2 to estimate the distance between the 
solution (y^,Z^) of BSDE(iV) (3.2) and that of the MFBSDE (3.3) we 
get, with the notation A = {X, Y, Z) and A^ = {X^ ,Y^ , Z^), that for any 
(5 E (0, 1) there are some 7 > and C > only depending on 5 (and, hence, 
in particular independent of N) such that 
rT 



E 



a"/*(\Y^ 



V|2 I I7JV 



Zty)dt 



<CE 



olT 



N 



N 



Y^{e^mx^,QHx^))-E[^{x,XT%^^, 



k=l 



+ 6E 



alt 



1 ^ 
i^(e'=<7)(i,Af,G^XAf)) 



N 



k=l 



E[g{t,u,At)]\^^At 



dt 



Hence, taking into account that g(u;,t,-,-) and $(a;,-,-) are Lipschitz with 
a Lipschitz constant L > which does not depend on {uj,t), we deduce that 

rT 



E 



5"*(ln 



N 



^/\2 ,\yN 



z.ndt 



<Rn + SCL^E[e^' \Xl^ - Xt^] 
+ 8L'6E I e^*(|A:f 



^t\ + K t 



Yt\' + \Z, 



N 



Zt\^)dt 
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1 ^ 

-Y.i®'<i>){XT,eHxT))-EMx,XT)]\ 



X=Xn 



fc=l 



=7* 







TV 



N 



Y^{@^g)(tAuQH^t))-E[g{t,n,At)]\^._ 



k=l 



dt 



Consequently, for 5 := (16L^) ^, 
rT 



2 



Jo 



<Rn + 8CL^E[e^' \X^ - Xxr] + -E 



re^'\X^-Xtfdt 
Jo 



E 



asiV^+cx). 



From Step 1 and (C3) it then follows that 



For proving that {Y ,Z ) converges to {Y,Z) also in B"^ we have still to 
show the convergence of Y^ to Y in 5p([0,T]). To this end we apply Ito's 
formula to jl^ — Yt\^ and take then the conditional expectation. Thus, be- 
cause of the boundedness of g, we have: 

\Yf -Yt\'^ + E 



J \Z^ - Zs\^ ds\J^t 



E 



1 



N 



N 
+ 2E 



rt 



Y,{B^^){X^,QHx^))-E[^x,XT%^Xr 

fc=i 

x|^E(0'5)(^>Af,e'=(Af))-i5;[5(.,u,A,)]|„=^^|(i.|.Fi 



fc=l 



<E 



1 ^ 
-Y.iQ^^)i^T,^\X^))-E[^{x,XT%._ 



:X7 



fc=l 



Tt 



+ CE 
Therefore, 
E 



f \Yf -Ys\ds\J^t 
Jo 



t€[0,T]. 



rN 



sup ly/^-Ftl' 

iG[0,T] 
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<AE 



1 ^ 



N 



k=l 



4n 



+ ^C^TE 



IK 



N 



YA'^ds 



and the announced convergence follows now from the preceding result and 
from Step 1. D 

3.2. Convergence speed of the approximation of the MFBSDE. In The- 
orem 3.1 we have seen that under the assumptions (C1)-(C3) the solution 
{Y^ , Z^) of BSDE(A^) (3.2) converges toward the unique solution {Y, Z) of 
our MFBSDE (3.3). The objective of this section is to study the speed 
of this convergence in the special case where the sequence {X^)i<[>i is 
an approximation of a forward SDE of McKean-Vlasov type. For this let 
6:0 X [0,r] X M'" X M'" ^ M™ and a:0 x [0,r] x M"^ x M™ ^ W^""^ be 
bounded measurable functions which are supposed to be Lipschitz in {x,x') 
with a Lipschitz constant C, that is, P-a.s., 

(D) \h{t, X, x') - b{t, X, x')\ + \a{t, x, x') - a{t, x, x')\ < C{\x -x\ + \x' -x'\), 
for all t€[0,T],{x,x'),{x,x')£R'^ xR"". 

We consider the following forward equation of McKean-Vlasov type 



(3.4) 



dXt = E[b{t,x,Xt%=x^dt + E[a{t,x,Xt)]i^=XtdWt, 

Xo = xeW, te[o,T], 

and we approximate this equation by forward equations in the same spirit 
as we have approximated our MFBSDE (3.3) by BSDE(iV) (3.2), A^ > 1. 
More precisely, we consider as approximating SDE(A^): 

1 ^ 



dXt"" = T7 Y.(®'b){t,x,'',e'{xn)dt 



N 



k=l 



N 



(3.5) 



X, 



N 



+ 4E(0M(t,^f,e^(xf))dw^„ 
fc=i 

X, te[o,T]. 



For the above forward equations we can state the following result. 

Proposition 3.2. Under the above standard assumptions on the coef- 
ficients b and a we have: 

(i) The forward SDE of McKean-Vlasov type possesses a unique solu- 
tion X £ 5p([0, T];M™). Moreover, X is adapted with respect to the filtration 
F generated by the Brownian motion W . 



16 



BUCKDAHN, DJEHICHE, LI AND PENG 



(ii) For all N > 1, the forward equation SDE(N ) (3.5) admits a unique 
solution X^ G cS|([0,r];M"). 
(iii) 



E 



sup \X^ 

.te[o,T] 



TV 



Xt\ 



< 



N 



for all N>1, 



where C is a real constant which only depends on the bounds and the Lips- 
chitz constants of the coefficients b and a . 

Proof. By using the properties of the operators Q'^^k > 1, the results 
(i) and (ii) can be obtained by easy standard estimates for SDEs. For this we 
see that the coefficients jiJ2k=ii®'^b) {t,x,x') and ;^I]^=i(0'^o")(t,x,x') of 
SDE(A^) (3.5) are bounded, F-progressively measurable for all {x,x'), and 
Lipschitz in {x,x'), uniformly with respect to {uj,t) £Q x [0,T]. Moreover, 
the bound and the Lipschitz constant do not depend on N. The coeffi- 
cients of the SDE (3.4) of McKean-Vlasov type, too, are bounded and Lips- 
chitz, uniformly with respect to t G [0,T]. From the fact that the coefficients 
E[b{t,x,Xt)] and E[a{t,x,Xt)] are deterministic it follows easily that the 
unique solution of the McKean-Vlasov type equation is F-progressively 
measurable. 

For the proof of statement (iii) of the proposition we notice that, thanks to 
the fact that the processes Q^{a{-,x,X)) = {{Q^a){t,x,Q^{Xt)))t(z[o^T],k> 
1, are mutually independent and identically distributed, of the same law 
as the process a{-,x,X) = io'{t,x,Xt))t£[o,T]^ ^^'^ independent of W (and, 
hence, also of X), we get 



E 



1 ^ 

-Y,{e'a){t,Xt,e\Xt)) - E[a{t,x,Xt)]^,=x, 



N 



k=l 



PxAdx) 



I ^ ^E(eM(t,^,0'(^t))-^k(i,^,^t 

E / y^4^T.(®'a,,){t,x,e\Xt))]pxM 

^ Y. [ y^ria,,,{t,x,Xt))Px,{dx)<^ 



l<i<Tn,l<j< 



l<i<m,l<j<d- 

for some constant C which does not depend on A'^ > 1. Similarly, we have 



E 



1 ^ 

-J2{&%){t,Xt,e'^{Xt)) - E[b{t,x,Xt)] 



N 



\x=Xt 



fc=i 



< 
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Using the fact that j^Ek=i{&''<^){tr,-) and jfEk=ii^''b){t,-,-) are Lips- 
chitz, uniformly with respect to t G [OjT], and we obtain now statement (iii) 
by an SDE standard estimate. D 

The above estimate of the convergence speed can be extended from the 
forward equations to the associated BSDEs. Indeed, we have the following. 



Theorem 3.2. We assume (C1)-(C2) and take instead of (C3) the as- 
sumption that X^ GcS2([0,r];M'^) is the unique solution of SDE(N) (3.5), 
N>1. Let {Y,Z) denote the unique solution of MFBSDE (3.3) driven by 
the forward SDE (3.4) of McKean-Vlasov type with solution X, and let 
(y^, Z^) he the unique solution ofBSDE(N) (3.2) driven by SDE(N) (3.5). 
Then, for some constant C which depends only on the bounds and the Lip- 
schitz constants of b, a and f , 



E 



sup \Yi 
.te[o,T] 



N 



Yt? + 



\Z-. 



N 



Zt\ 



(J 

< — for all N>1. 

~ N 



Proof. We first notice that, since due to Proposition 3.2 the process 
X is F -progressively measurable, the unique solution {Y, Z) of MFBSDE 
(3.3) driven by the process X must also be F -progressively measurable. 
Consequently, the triplet A = (X, Y", Z) is independent of all the processes 
G'=(/(-,u,A)) = ((e*^/)(i,u,G^(At)))i6[o,T], k>l, which allows to use the 
argument developed in the proof of Proposition 3.2 and to show that, for 
some constant C which only depends on the bound of /, 



E 



1 ^ 
-5:(eV)(i,Ai,e'=(A,))-i?[/(t,u,A,)] 



|u=At 



fc=l 



- N 

for ah iV > 1. 



Analogously, we get, again for some constant C which only depends on the 
bound of <1>, 



N 



Y,{Q''^){XT,<d\XT))-E[^{x,XT)]\.=X, 



k=l 



2n 



< 



N 



for all N>1. 



These estimates together with Proposition 3.2 and BSDE standard estimates 
yield the wished speed of the convergence of {Y ,Z ) to (Y,Z). The proof 
is complete. D 



4. A central limit theorem for MFBSDEs. We have seen in the preced- 
ing Section 3.1 that, if the sequence X ,N > 1, is defined by the forward 
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equations SDE(A^) (3.5), the speed of the convergence of {X^ ,Y^ ,Z^) to 
{X,Y, Z) is of order 1/^/N, that is, for some real C, 



E 







/•^ .r 


- 


sup(|Xi^- 


- ^t\ +\yt - 


-yt?)+ i^f- 


-Zt\^dt 


te[o,T] 




Jo 


- 



<— forAf>l. 

- N 



Consequently, the sequence of processes y/N{X'^ — X,Y'^ — Y,Z^ — Z), 
N>1, is bounded in cS|([0,T];M'") x S^{[0,T]) x £|([0,r];M^) and the 
question about the limit behavior of this sequence arises. The study of this 
question is the objective of this section. For this end we will begin to in- 
vestigate in a first section the limit behavior of the sequence (viV(X^ — 
X))n>i C Sp{[0,T];MJ^). This discussion which involves, for its most part, 
some recall from known facts prepares the study of the limit behavior of the 
triplet VN{X^ - X,Y^ - Y,Z^ -Z), N> 1, which wih be the object of 
the second section. 

For the sake of simplicity of the notation but without restricting the 
generality of our method we suppose in what follows that the coefficients 
6,0",/ and $ do not depend on {uj,t). Moreover, we will assume that the 
functions are continuously differentiable. 

4.1. Limit behavior of y/N{X^ -X). Let 6 : M"^ x M'^ ^ R'^ and a:R'^x 
M — > M be two bounded, continuously differentiable functions with bounded 
first-order derivatives, and let xq be an arbitrarily fixed element of M'^. As 
in the preceding section, but now under the additional assumption that the 
coefficients do not depend on {co,t), we denote by X'^ G 5p([0, T];]R'') the 
unique solution of SDE(A^) 

rt 1 N 



(4.1) 



^0 ^ k=l 

rt 1 N 



fc=i 

-.9 



and X G 5pvi.'([0,T];M ) is the unique solution of McKean-Vlasov equation 

Xt = xo+ f E[bix,Xs)]i^=Xsds+ f E[a{x,Xs)]i^=XsdWs, 

(4.2) ^° ^° 

s€[0,T]. 

Then we have the limit behavior of (//V(X^ - X))n>i C cS|,([0,r];M'^): 

Theorem 4.1. iei^ = (?«,?(')) ={((d'"^(^))i<.<d,(d'"'''^(x))i<„<d), 
{t,x) G [0,r] X W^} be a d + dx d-dimensional continuous zero-mean Gaus- 



MFBSDES: A LIMIT APPROACH 19 

sian field which is independent of the Brownian motion W and has the co- 
variance functions 

E[^l'''\x)cl?''^'\x')]=cov{b\x,Xt),a'^'{x',Xt>)), 

E[d'''''\x)cl?''^'\x')]=cov{a^'^{x,Xt),a''\x',X,)) 

for aU{t,x),{t',x')e [0,T] xR'^,1 <i,j,k,£ < d, 
and let X 6 6':2.Q0,T];]R ) be the unique solution of the forward equation 

Xt= f6'Hx,)ds+ fif\Xs)dWs 

JO JO 

(4.3) + [ {E[{VMx.Xs)]a,=xXs + E[{V.A){x,X,)X,]^=x^ds 

Jo 

+ / {E[{V,,a){x,Xs)]^=x.Xs + E[{V,,a){x,Xs)Xsh=xJdWs, 
Jo 



t G [0,T], where F is the filtration F augmented by a{£,tix),{t,x) G [0,r] x 
W^}. Then the sequence {\fN {X'^ -X))m>i converges in law over C{[0,T];W^) 
to the process X. 

Remark 4.1. The continuity of the above introduced two-dimensional 
zero-mean Gaussian process ^ is a direct consequence of Kolmogorov's Con- 
tinuity Criterion for multi-parameter processes. 

Indeed, a standard argument for mean-zero Gaussian random variables 
and a standard SDE estimate shows that, for all m > 1 and for some generic 
constant Cm which can change from line to line, 

E[Mx)-iAx')\''^] 

<CmJZE[\i^^''\x)-i\}^\x')\'-] 

1=1 
+ CmJ:E[\e^'^^\x)-if''^^\x'f-] 



<CmY.iE[\e'\^)-^'v'\^')\"]T 
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d 
= CmY.{\^T{b\x,Xt) -b\x' ,Xt,))r 



i=l 

d 



<C„(|t-t'r + |x-x'p'") for all (t,rE), (t',x') G [0,r] x M. 

The proof of Theorem 4.1 will be split into a sequel of statements whose 
objective is the study of the limit behavior of the process \/N{X^ — X) 
as A^ tends to infinity. The process y/N{X — X) can be described by the 
following SDE: 

Vn{x^ -Xt) 

= j\es'''{X^) + VN{E[b{x,X^)]^,^^^ -E[b{x,X,)]\,^xJ}ds 

+ j\es'''{X^) + ^{E[a{x,X^)]\,^x^-E[a{x,X,)]\,^xJ}dWs, 

t£[0,T], 



where 



TV 



e^^(x) = J2ib{x,&Hxn) - E[b{x,Xl')]), 



Xk=i 



N 



1 

i}''{x) = -^Y.(a{x,^\X^))-E[a{x,Xf)]) 



k=l 

{t,x)e [0,T] xM'^. 

Obviously, the random fields f'^ = {C'^ix), {t, x) G [0, T] x M'^}, i = 1, 2, 
are independent of the Brownian motion W and they have their paths in 
C([0,r] X W^;W^) and C{[0,T] x M'^;M'^^'^), respectively Moreover, we can 
characterize their limit behavior as follows: 

Proposition 4.1. The sequence of the laws of the stochastic fields ^ = 
{^^'^,^^'^),N > 1, converges weakly on C([0,T] x R'^;]R'^ x M'^^'^) toward 
the law of the continuous zero-mean Gaussian (d + l)-parameter process 
C = {{S,^^,^^'^')} introduced in Theorem 4-1- 

The proof of the weak convergence of the laws P o [,^"]~^, A^ > 1, on 
C([0,T] X M'^;E'' X M'^^'^)) is split into two lemmas: while the first lemma 
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establishes the tightness of these laws on C([0, T] x M*^; R'^ x M'^^'^) the second 
lemma will study the convergence of their finite-dimensional marginal laws. 
The both results together, the tightness of the laws and the convergence of 
the finite-dimensional marginal laws imply the weak convergence. 



Lemma 4.1. The sequence of the laws of the stochastic processes ^ 



N 



(C^'^,^'^), N>1, is tight on C{[0,T] x 



ixd 



) . In fact, we have 



for all m>l the existence of some constant Cm such that: 

(i) ^[|Cf (x)|2-] < Cm., for all {t,x) G [0,T] x M'^, A^ > 1; 
(ii) E[\^^{x) - ei/(x')|2-] < Cm{\t - t'r + \x- xf-), for all {t,x), {f , 
x')e[0,T] xR'^,A^>l. 

Moreover, for all m>l there exists some constant Cm such that 



(hi) 



E 



sup 



1 + Ixl 



2m. 



<c„ 



for all t£ [0,r],iV> 1. 



Proof. We begin with the proof of statement (i). Putting 



-i,N 



^N\ 



b"''{t,x):=b\x,Xn-E[b\x,Xn] 



and 



a''^^^{t,x) := a'^^{x,X^) - E[a''^{x,X-, 



N\ 



l<i, j <d, N>1, 



we have the existence of some generic constant C^^m such that, for all (i, x) G 
[0,r] xM*^ andfor aU 7V>1, 



E[\^f{x)\'"^]<Cm4Y.E 



\i=l 



N 

Y.QHV^''{t,x)) 



k=l 



2m- 



+ Y.E 
«J=1 



N 



)=Y^Q\a-'-^{t,x)) 



k=l 



2m 



On the other hand, taking into account the independence of the zero- mean 
random variables Q^{b ' (t, x)), 1 < k < N and the boundedness of the func- 
tion b we get, with the notation Tm,N = {(^ii • • • i ^2m) G {1, . . . , A^}^'": for all 
«(1 < ^ < 2m) there is a j G {1, . . . , A^} \ {i} s.t. ki = kj}, 



E 



1 



N 



j:eHV'\t,x)) 



fc=l 

1 



2m- 



N 



m' 






■2m 

ne^-^(F'^(i,x)) 

■i=i 
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"2m 



N^' 



E 



E 



{k-l,...,k2m)&Tm.,N 



X{e'K^'\t,x)) 



■i=i 



2m 



< 



N^: 



y: n^[i0'^(^''' (*'-))! 



|2m]l/(2m) 



(fcl,...,fc2m)Grm,]vi = l 



<^card(r„,Ar)^[|6*' (t,x) 



2mi 



\,N 



for some C^, C^ G IR+, which are independent ofA^>l, 1 < i < d and 
(t,x) e [0,r] X M'^. Similarly, we see that, for all iV > 1, l<i,j <d and 

(t,x)G[0,r] xM'^, 

2mn 



^ 



AT 






<c„ 



Consequently, for some Cm G IK+ neither depending on A^ > 1 nor on (t, x) G 
[0,r] X M'', £J[|^j^(x)p™] < Cm, and the same argument, but now applied to 
V.6(x,x'), V,a(x,x'), also yields E[\V ^^^ {x)\^^] < Cm- 

For proving statement (iii) of the lemma we observe that, for all m > 
{d+ l)/2, due to Morrey's inequality and the estimates obtained above, we 
have the existence of some generic constant Cm depending on m and d such 
that, with the notation ^f (x) = (1 + |x|)-^^f (x) and 7 = 1 - 5^, 



E 



sup 



l + lxl 



2m- 



<E 



( sup |fi(x)| + sup 



Ux)-Ux')\ 



< CmE 



<a 



<C« 



X — x'l"^ 

W/„M2m , |V7 ^Nf„\i2m\ 



2m 



(|Cf(x)|^™ + |V,Cf(x)|''")dx 



iEm\x)n+E[\v,^^ixr"^]) 



dx 



(1 + kl 



\2m 



for all t G [0,T],N > 1. To get (ii) we follow the argument given for the proof 
of (i). Combining it with a standard SDE estimate we get, for a generic 
constant Cm and ah (i,x), {t',x') G [0,r] x M.'^ and A^ > 1, 

E[\^^{x)-^ffix')\'n 

/ d 



i,N , 



<Cm\Y.E{\\}''\t, 



r^(t',x')pn 



^i=l 
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i,7 = l / 



<Cm{E[\b{x,X(')-b{x',X(/)\'^]+E[\a{x,X[')-a{x',X^)f^]) 

<Cm{\t-t'r + \x-x'\^"'). 

Finally, to complete the proof it only remains to observe that due to Kol- 
mogorov's weak compactness criterion for multi-parameter processes the es- 
timates (i) and (ii) imply the tightness of the sequence of laws of ^ , A^ > 1, 

on C([0,r] xM'^;M^xM'^>''^). D 

For proving Proposition 4.1 we have still to show the following lemma. 

Lemma 4.2. The finite- dimensional laws of the stochastic fields ^ = 
(^"^' , ^^' ), A^ > 1, converge weakly to the corresponding finite- dimensional 
laws of the continuous zero-mean Gaussian {d -\- l)-parameter process ^ = 
{{Ct^\x),Ct^\x)),{t,x) G [0,r] X W^} introduced in Theorem 4.1. 

Proof. We decompose the random field S,^ = (C^'^,C^'^). For this we 
first consider the component ^^' and we represent it as the sum of the 
random fields ^^'"'^ and ^^' '^, where 

1 ^ 
^^''^'^(x) = ^(e'=(cT(x,Xi)) - E[a{x,Xt)]), {t,x) e [0,r] X M^ 



1 ^ 

for {t, x) G [0, T] X R'^. Since the bounded {d+ l)-parameter fields e^{a{-,X)) -- 
{a{x,Q^{Xt)),it,x) G [0,r] X M.'^}, k>0, are i.i.d., it follows directly from 
the central limit theorem that the finite-dimensional laws of the random field 
^2,N,i converge weakly to the corresponding finite-dimensional laws of the 
zero-mean Gaussian {d -\- l)-parameter process ^'^' whose covariance func- 
tion coincides with that of the field ct{-,X) = {a{x,Xt), {t,x) G [0,r] x M"^} 
[recall the definition of .^ = {C ,C ) given in Theorem 4.1]. On the other 
hand, since due to Proposition 3.2, 

1 ^ 

= -y E 



N 




'y]{e\a'^\x,Xf)-a'^\x,Xt)) 


k=i 




2- 


-E[a''^{x,X[')-a''^ix,Xt)]) 
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c 

<— for alliV>l,(t,a;)G [0,T] xM'^, 

it follows that also the finite-dimensional laws of the random field ^^' = 
^2,N,i _|_ ^2,^,2 converge weakly to the corresponding finite-dimensional laws 
of the zero-mean Gaussian field ^^"^^ . 

By applying now the argument developed above to the couple ^'^ = 
(C^ iC^ ) we can complete the proof. D 

As already mentioned before, Proposition 4.1 follows directly from the 
Lemmas 4.1 and 4.2. Proposition 4.1 again allows to apply Skorohod's rep- 
resentation theorem. Taking into account that the random fields ^^, A^ > 1, 
are all independent of the Brownian motion W, we can conclude from Skoro- 
hod's representation theorem that, on an appropriate complete probability 
space {ft' ,J-'' ,P') there exist {d + dx d) -dimensional {d+ l)-parameter pro- 
cesses C'^ = {^^^(x),(t,x) E [0,r] X M'^} {N > 1) and ^' = U^(x),(t,x) G 
[0,T] X M*^}, as well as a d-dimensional Brownian motion W' = (W/7)te[o,T]) 
such that: 

(i) Pl,=P^,P^,^=P^M,N>l; 

(ii) ^'^ = (^'i'^,^'2,7V) — ^ ^/ ^ (?'^?'^), uniformly on compacts, P'-a.s.; 
(hi) W is independent of ^' and (,'^ , for ah iV > 1. 

For this new probability space we introduce the filtration F' = (JT^ = J^j^ V 
•n)ie[o,T], where ^^ = aU^(x),Cf (x), (s,x) e [0,r] x R^iV > 1} V AA, and 
we observe that W is an F'-Brownian motion. 

Given £^'^ , ^' and the Brownian motion W' we now redefine the processes 
X^ and X on our new probability space. For this we recall that the process 
X^ is defined as unique solution of (4.1). In virtue of the definition of 
^^ = (C^'^,'^^'^) we can rewrite (4.1) as follows: 

Xf = ^o + J^ ^-^es'^'iX^) + E[b{x, Af )]|,=;,^.) ds 

+ J^(^^es'''{X^) + E[a{x,X^)]^,^,,MyWs, te[0,T]. 

Taking into account that, P-a.s., 

\0^ix)-C'^{x')\ <2VnL\x-x'\ for ahsG [0,T],x,x' eM.'^,i = 1,2, 

where L denotes the Lipschitz constant of b and of a, the fact that ^^ = 
(^^' , ^^' ) and ^ = (C'^ )C'^ ) obey the same law has as consequence 
that also C^^' (•) and £,'^' {■) are Lipschitz, with the same Lipschitz constant 
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as S,l'^ and ^g''^ ■ Therefore, the equation 
Xr = xo + J^ (-^e^^(Xf ) + E'[b{x, Xf )]|,=^,.) ds 



admits a unique solution process X'^ S 5p/([0,r];M'^). In the same spirit in 
which we have associated with X the process X we define the analogue 
to X on the probability space {n',J='',P') by letting X' £ 5|,([0,r];M'^) be 
the unique solution of the SDE 

(4.4) ^° ^' 

tG[o,r]. 

We remark that the solution X' is adapted to the filtration generated by 
the Brownian motion W' . Moreover, the fact that the laws of (^ , W) and 
(^ , VF) coincide implies that also 

p' o (X', x'^, e'^, w'Y^ = Po{X, x^,c^, W)~\ N>1. 

Consequently, Theorem 4.1 follows immediately from the following proposi- 
tion: 



Proposition 4.2. Let X G 5|,([0,r];M'^) be the unique solution of the 
SDE 

Xt = j\i'hX's) + E'[V.Mx,X',)]\,^x'Xs 

+ E'[V^,,b{x,X',)Xs\\.=xOds 
(4.5) 

+ ^ {i'!{X',) + E'[VM^,K)]\.=X'Xs 

+ E'[V,.a{x,K)Xs\\^=xOdK, 
t€ [0,r]. Then 

E' 



.se[o,t] 



sup \y/N{X'f -X',)-Xsi 



as N ^ +00. 



Proof. We begin by showing 

Step 1. The existence and uniqueness of the solution X £ 5f'([0,T]; 
of the above SDE (4.5). 
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For this we remark that, thanks to the independence of ^' = {^'^,C''^) and 
W' we have also that of ^' = (^'\C'^) and X' , so that 



<IUt:lif vl\\2mi 



E'mix. 



E'[\C{x)r]P'x'{dx)<C„,, sG[0,T], m>l 



(cf. Lemma 4.1). Since the other coefficients of the above hnear SDE are 
bounded we have the existence and the uniqueness of the solution X in 
5F'([0,r];M ); moreover, X is adapted with respect to the filtration gener- 
ated by W and ii^'[supjgro ^-i \Xtf] < oo, for all p > 1. 

Step 2. (er(X^^))te[o,T] -^ (C{(^t))i6[o,T] m L\[Q,T] x n',dtdP'), as 
Indeed, for all ?n > 1, we have 

E'[\C{x'j')n=E[\c^{x^)r] 



< E 



sup 

■x£Rd 



l + lxl 



2m-, ^ 1/2 



'(i?[(i+ixfm)V2. 



Thus, using Lemma 4.1(iii) and the fact that the coefficients of (4.1) are 
bounded, uniformly with respect to A^ > 1, we can conclude that, for some 
constant Cm, 

E'[\C{Xf)n < Cm for all s G [0,r],iV > 1. 

On the other hand, recalling that ^'^ = (C'^'^,C'^'^) ^ ^' = (C'^C'^) uni- 
formly on compacts of [0, T] x W^, P'-a.s., and 



E' 



sup \X'f-Xi\' 
.se[o,T] 



E 



sup \X^-Xsf 

SG[0,T] 



c 



N>1 (cf. Proposition 3.2), 

we see that sup^gp^T] ICJ^ i^'J^ ) ~ ^si-^s)\ ~^ 0' i^ probability. Combining this 
with the uniform square integrability of {^'^^ {Xg^),s G [0,T],N > 1}, proved 
above, yields the convergence of^'^{X'^) to ^'{X') in L'^i[0,T] x 9.',dtdP'). 



Step 3. Let us now decompose \/N{X[" — X[) — Xt as follows: 



iv(xf - x[) -Xt = P/' + P/' + It : i e [0, T] 



2,N I j3,N 



where 



il'^= f\c'^^^^ixf)-eHxi))ds+ f{C'^{x'^)-i?{K))dw',, 

Jo Jo 

It"" = j\VN{E'[b{x,Xf)]\,^^,. - E'[bix,Xi)]\,=xO 



t3,N 
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-iE'[VMx,X',)]\,=X'Xs + E'[VMx,X',)Xs]\,=x'J}dW',. 

From the Step 2 we know already that i?'[supjgrQ ^i |/j ' p] ^ 0, as A^ — > +oo. 
Let us now esthnate £^'[sup^gro ti l-^t ' P]- ^^^ this we remark that, thanks to 
the continuous differentiabihty of the function a, with the notation X'^ '^ := 
X^ + 7(Xf -X^)weget 

rt ( /•! 



r3,Af 



L^O 



£;'[V.a(xi,Xf '^) - V,.a(x2,X;)], _^,iv,. ^_^, d^ x X 



ixi=xr'^,^2=x^ 



JO 



V^'a(x2,X^)}X,]|^^^^,^,^^^^^^, dTJ dP^i 



+ y^ {E'[V.cT(x,Xf '^)]^^^,^„(^/iV(Xf - X',) - Xs) 

+ E'[V^,a{x, X'f'^^iVNiX'f - Xi) - Xs)l^x'f,.}dWi, 

t£[0,T], 

and taking into account that the first-order derivatives of a are bounded 
and continuous, it follows from the convergence of X"^'"^ — > X' (N ^ +00) 
in 5p/([0,T];]R'^), for all 7 G [0, 1], that for some sequence < pat \ 0, 



E' 



sup |/; 

se[o,t] 



3,Ni2 



<PN + C E'llVNiX'/" -Xi)-Xsr]ds 



t£[0,T], N>1. 
The same argument yields the above estimate also for /^'^. Consequently, 

E' sup \VN{X'f - Xi) - Xs\^ 
■se[o,t] 



<E' 



sup I If 

i6[0,r] 



l,Af|2 



+ 2pN + 2C E'[\VN{X'f -X'J-XsHds, 



te[0,T], N>1, 



and from Gronwall's lemma we obtain 



E' 



sup \VN{X'f - X',) - Xs 
se[o,t] 



as A^ ^ +00. 
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The proof is complete. D 

It remains to give the proof of Theorem 4.1. 

Proof of Theorem 4.1. In analogy to the existence and uniqueness 
of the solution X € 5|,([0,r];M'^) we can prove that of X G 5|,([0,r];M'^). 
Moreover, since W is independent of ^ and ^ ,N > 1, and, on the other 
hand, also W' is independent of ^' and ^ ,N > 1, it follows that: 

(i) Po{W,^^)-^ = P'o{W',^'^)-\ for aniV>l; 

(ii) Po {W,i)-'^ = P'o {W',^')-^. But X'^ is the unique solution of an 
SDE governed by the W and £,'^ , and X^ is the unique solution of the 
same equation but driven by W and ^^ instead of W and ^'^. On the 
other hand, X' is the unique solution of (4.2) with W' at the place of W. 
Thus, 

(iii) Po{W,^^ ,X,X^)-^ = P'o{W',^'^,X',X'^)-\ioi all N >1. Con- 
cerning the process X, it is the unique solution of the same equation as X, 
but with ^ and W replaced by ^' and W , respectively. This, together with 
(ii), yields _ 

(iv) Po {W,C,X,Xy^ = P' o {W',C,X',X)-^. From (iii) and (iv) we 
see, in particular, that P^(^xn~x) = ^[/NiX'^-X'Y ^ - ^' ^""^ ^'^^^'x 
For completing the proof it suffices now to remark that the convergence 
of ^/N{X'^ - X') to X in cS|,([0,r];M'^) implies that of their laws on 
C([0,T];M^). D 

4.2. Limit behavior of VN{X^ - X,Y^ - Y,Z^ - Z). As in the pre- 
ceding section we suppose that 6 : M*^ x M"' ^ M'' and cr : M"^ x M*^ ^ M*^^ '^ 
are bounded, continuously differentiable functions with bounded first-order 
derivatives, and xq is an arbitrarily fixed element of M . With the forward 
equation 

rt 1 N 

X^ = xo+^Y.b{X^^@\X^))ds 
•^0 ^ k=l 

^0 ^ k=l 

we associate the BSDE 

N „T -, N 



^t"" = ^ E ^(^T , Q'ix^)) + / ^ E /(Af ' 0'(Af )) ds 

k=^ k.=^ 

(4.7) 



T 
t 



ZfdWs, 
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where 

A^ = (X^,y^,Z^)G5^:=cS|([0,T];M'^)xcS|([0,r];M)x4([0,r];M'^) 

is composed of the unique solution X of the forward SDE and the unique 
solution (y^, Z^) of the BSDE. The functions / : (M"^ x M x M'^)^ -^ R and 
$ : R*^ X W^ -^ M are supposed to be bounded and continuously differentiable 
with bounded first-order derivatives. Moreover, our approach imposes the 
following additional hypothesis for the function /: 

(E) f{{x,y,z),{x',y',z')) = f{{x,y,z),{x',y')), for all {x,y,z), {x'y',z') 
e M"^ X M X M"^. 



We know already that 



E 



sup {\x^ -x,\' + \Y,^ - y,|2) + r |zf - z, 



<§, N>,. 



for some constant C, and we are interested in the description of the limit 
behavior of the sequence \/N{X — X, Y — Y, Z — Z). A crucial role in 
these studies will be played by the following uniform estimate for Z^^ ,t G 

[o,r],A^>i. 

Lemma 4.3. Under the above assumptions on the coefficients there exists 
some real constant C such that, dtdP-a.e., 

\Zt\ < C and |Zf | < C for all N>1. 

Proof. We remark that A^ = {X^ ,Y^ , Z^) £ B^ is the unique solu- 
tion of the system 

X^ = xo+ fb^{X^)ds+ fa^{X^)dW,, 
Jo Jo 

Yf = <^^{X^)+rf^{A^)ds- r Z^dW,, tG[0,r], 
Jt Jt 



where the coefficients 

k=i k=i 

, Af ^ N 

fe=l fc=l 

t G [0, T], {x, y, z) G M'^ X M X W^, are bounded random fields such that: 

(i) (T^(a;),6^(2;),$^(x) and f^'{x,y,z) are .7^o-™easurable, for all (i,x, 
y,2;), and, hence, independent of the driving Brownian motion W] 
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(ii) a^{-),bf{-),^{-) and /t (•) are continuously differentiable and their 
derivatives are bounded by some constant which does neither depend on 
t € [0, T] nor on w G il. 

Let Dt = {Dj, . . . ,Df),t G [0,r], denote the Mahiavin derivative with re- 
spect to the Brownian motion W = {W^, . . . , W^). For a smooth functional 
F G 5 of the form F = Cip{W^^ ,..., W/;), with n > 1, ii, . . . , t„ G [0, T], 1 < 
ii,...,in<n,(pe C^iR'"'"^) and C G L°°(fi, J^o, P), the MaUiavin derivative 
DfF is defined by 

n 

The operator D:S C L'^{n,J=',P) -^ L'^{n x [0,T],dPdt;R'^) is closable; its 
closure is denoted by (D,D"^'^). 

It is well known that, for the coefficients b ,a , f and <^ with the 
above properties, X^ ,Y^ and Z^ belong to L^{[0,T];]n)'^'^), Xf^ ,Yf G B^'^, 
for ah tG [0,r], and 

J s J s 

DiYf = DiX^V,<^''{X^)+rDlA^Vxfr''{A^)dr- j^ D\Z^ dWr, 

Jt Jt 

0<s<t<T,l<i<d, where <7^'*'' denotes the ith row of the matrix a^ . 
The above result is standard in the Malliavin calculus, the interested reader 
is referred, for example, to the book of Nualart [9]. From standard SDE 
estimates we then get 



E 



sup iDsYt^l^lJ's 



<C, P-a.s., for alltG [0,T],iV>l. 



Ue[s,T] 
On the other hand, differentiating the equation 

Jo Jo 

in Malliavin's sense we obtain, for < s <t <T and 1 <i < d, 

DlYf = Zf^'- f'DlA^f,''{A^)dr+ f DlZ^ dW^. 

J s J s 

Consequently, letting s \ t we have DgY/^ ^ Z^ in L'^{Q,J^,P), ds-a.e. 
From there and the above estimate for D\Yi we deduce that 1^^ I — ^) 
dsdP-a.e., for all A'^ > 1. Finally, since this estimate is uniform with respect 
to A^, it holds also true for the limit Z of the sequence Z^ , N >1. The proof 
is complete. D 
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From the above lemma and the assumption of boundedness of the coeffi- 
cients ^ and / of the BSDEs we get the following by a standard estimate: 

Corollary 4.1. For all m>l there is some constant Cm such that, 
for all t, t' G [0,r] and all N>1, 

E[\Yf - y/pn < Cm\t - t'r. 

The same estimate holds true for the limit Y of the sequence Y ,N > 1. 

Recall that the limit X of the sequence X^ ,N> 1 , is the unique solution 
of the forward equation 

ft rt 



(4.8) 



Xt = xo + j E[h{x,X,%=Xsds + j E[a{x,Xs%=XsdWs 



tG[0,T], 

while the limit of the sequence {Y^ ,Z^) is given by the unique solution 
(y, Z) of the BSDE 

Yt = Emx,XT)]^=XT+ I E[f{X,As)]x=A.ds- f Z^dWs, 
, , . Jt Jt 

(4.9) 

tG[0,T], 

where A = (X,y,Z)GS|. 

We have the following limit behavior of {y/N{X^ - X,Y^ - Y, Z^ - 
Z))n>i C S| . 

Theorem 4.2. Let C = (e«,^(2)^^(3)^^(4)) ^ mi'^'\x)),<,<d, 
(Cp'^')(x))i<,,,-<rf,e(3)(x), ef)(x,y,z)), {t,x,y,z) e [0,T] xM.'^ xRx R'^} be 
a {d + dx d + 2) -dimensional continuous zero-mean Gaussian field which is 
independent of the Brownian motion W and has the covariance function 



E 



i^^\x) 

Ui'\x,y,z)J 



(S) 



i'^^Xx') 

\Cif\x',y',z')J 



E 



I b{x,Xt)-E[b{x,Xt)\ \ I b{x',Xt>)-E[b{x',Xt>)]\ 







a{x,Xt)-E[a{x,Xt)] 
<^{x,Xt)-EMx,Xt)] 

V [/(A,t)] ; V 

where [f{\,t)]:=f{\,At)-E[f{X,At)], X = {x,y,z),X' = {x' ,y',z'), {t,x, 
y,z),{t',x',y',z')e [0,T] x M'^ x M x M'^. 



aix',Xt')-E[aix',Xt>)] 

<^{x',Xt)-EMx',Xt)] 

[/(A',tO] 
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Recall that, for a = (a*)i<i<m, b = (fe*)i<i<m. G I^"^; a®h = {a^V)i<i,j<m. £ 
M"^x";Cf (x) = (Cp'^')(x))i<,j<rf and a{x,x^={a^i{x,x'))i<,,,<d are re- 
garded as d'^ -dimensional vectors. By A= {X,Y,Z) G B^ we denote the 
unique solution of the system 

Xt= fci'\Xs)ds+ f\f\Xs)dWs 
Jo Jo 

(4.10) + f\E[{V.Mx,Xs)].=x^Xs + E[{V,,b){x,X,)Xs].=xJds 

Jo 

+ / {E[{V,a){x,Xs)]c,=x.Xs + E[{V,,a){x,Xs)Xsh=xJdWs, 
Jo 

Yt = {C^^\XT)+E[V^^X,XT)h=XrXT + E[V,Mx,XT)XTl=Xr} 



(4.11) * _ _ 

+ E[Vxf{X,As)]x=AAs + E[^yf{\As)As]x=Ajds 



J ZsdWs, t£[0,T]. 



Recall that F is the filtration F augmented by the a-filed generated by the 
random field S,. Then the sequence {y/N{X^ - X,Y^ - Y,Z^ - Z))n>i 
converges in law over C([0,r];M'^) x C([0,T];M) x L'^{[0,T];W^) to A = 
(X,Y,Z). 

Remark 4.2. As in the case of the only forward equation (cf. Remark 
4.1) we get the pathwise continuity of the above introduced zero- mean Gaus- 
sian field with the help of Kolmogorov's Continuity Criterion for multi- 
parameter processes. From Remark 4.1 we know already that for all m > 1 
there is some constant Cm such that 

mfHx) - #(x')i'"i < Cm{\t - t'r + k - xfn 

for all {t,x),{t' ,x') E [0,T] x M , i = 1,2. The same argument also shows 
that, again for some constant Cm, 

E[\^^^\x) - ^(3)(2;')|2™] < Cm\x - xf'" for ah x,x' G M.'^. 
Moreover, for some generic constant Cm we have 

i^[id'^(x,y,z)-ef(x',y',z')i'™] 

<Cm{EUi'\x,y,z)-^f{x',y',z'f]r 

= CmiE[yarifi{x,y,z),At)-fiix',y',z'),At>))]r 
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< C^{\x - xf ™ + \y- yf ™ + \z- zf ™ 
+ E[\Xt-Xt,\^'^ + \Yt-Yt:\^'^]) 

for all t,t' £ [0,T],{x,y,z),{x',y\z') G M'^ x M x M'^ [recall the assumption 
(E)]. From standard SDE estimates we have 

E[\Xt-Xt>\^^]<Cm.\t-t'r, t,t'e[0,T]. 

Since, on the other hand, due to Corollary 4.1, we also have 

E[\Yt-Yt>\'^"']<C^\t-t'r for ant,t'G [0,T], 

we can conclude that 

E[\^l'\x,y,z)-^i^\x',y',z')fn 

< Cm{\t - t'r + k - a;f "^ + \y- yf "" + \z- z'f"") 

for all t, t' £ [0, T] , {x, y, z), {x' , y', z')eM'^xRx R°'. 

The proof of Theorem 4.2 is split in a sequel of statements which translate 
the proof of Theorem 4.1 into the context of a system composed of a forward 
and a backward SDE. For this we note that we can characterize \/N{Y — 
Y, Z — Z) as unique solution of the following BSDE: 

y/N{Yl' -Yt) 

= e'''{X^) + VN{E\^{x,X^)]^^^.-E\^x,XT)].=Xr) 
rT 



+ 1 Uf'^(Af) + ViV(i?[/(A,Af)],=A^-i?[/(A,A,)],=Aj}d. 

-J ^{z^-Zs)dWs, te[o,T], 



where 



1 ^ 



e'^'ix) = -j= Y^{Hx,Q\X?)) - E[^{x,X^)]), 

(t,x)G[0,r]x 



^.=1 



1 ^ 
ef'^(A) = ^=^(/(A,e^(Af))-£;[/(A,Af)]), 

V ^^ fc=l 

(t,A)G [0,T] y.W^y.Ry.W^. 

Recall also the definition of the random fields ^*'^ = {^^' {x),{t,x) G 
[0,r] X M"^}, i = 1,2, which have been used for rewriting the (forward) 
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equation for X . We remark that the random field ^'' = {^l' (x), {t,x) G 
[0,r] X M'^}, i = 1,2, e^'^ = U3,iV(x),x G M'^} and ^^'^ = {^'^ {X),{t, \) G 
[0,T] X M'^ X M X M'^} are independent of the Brownian motion W and their 
paths are jointly continuous in all parameters. Moreover, we can characterize 
their limit behavior as follows: 

Proposition 4.3. The sequence of the laws of the stochastic fields ^ = 
(?^'^)C^'^)^^''^)^'^'^))-^ ^ 1) converges weakly on 

C := C([0,r] X R'^;R'^) x C([0,r] X IR'^;IR'^>''^) x C(M'^;M) 

xC([0,T] xM'^xIRxM'^;]R) 

to the law of the continuous zero-mean Gaussian field ^ = {^^^' , C , C ^\ ^ ) 
introduced in Theorem J^.2. 

The proof of the proposition uses the same approach as that developed 
for the proof of Proposition 4.1. In a first lemma we establish the tightness 
of the sequence of laws PtN ^ N > 1, on C, and a second lemma is devoted 
to the convergence of the finite-dimensional distributions of the sequence 
{^ , A'^ > 1}. Both statements together yield the weak convergence of the 
laws PfN, N >!, on C. 

Lemma 4.4. The sequence of the laws of the stochastic processes ^ = 
(^^i,N^^2,N^^3,N^^4,N^^ A^ > 1, is tight on C. In fact, in addition to the es- 
timates established in Lemma 4-1 for ^*' , A^ > l,i = 1,2, we have for all 
771 > 1 the existence of some constant Cm such that: 

(i) E[\f^^{x)\'"' + |ef'^(x,y,z)|2-] < Cm 

for all{t,x,y,z)e[0,T] x M"^ x E x E'^,iV > 1; 

E[\f'^{x)-f'^{x')r + \^P''{x,y,z)-^t'"'i^',y'^^')n 
< Cm.{\t - t'r + \x- xf" + \y- yf "^ + \z- zf"), 

for all{t,x,y,z),{t',x',y',z')£ [0,T] x M'' x M x M'',iV > 1. 

Moreover, for all m>l there exists of some constant Cm such that 



(iii) E 
for allte [0,T],N> 1 



sup 1 I I I I I I 



<c„ 



Proof. The proof of the estimates (i) and (iii) uses the same arguments 
as those developed in the proof of (i) and (iii) of Lemma 4.1, with the only 
difference that now Morrey's inequality has to be applied for m > d + 1 and 
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J = I — 2^ ■ ^s concerns the proof of estimate (ii) we follow the argument 
developed for the proof of (ii) of Lemma 4.1. So we get, for a generic constant 

Cm and all (t,x,y,z), (t',x',y',z') S [0,r] x M'^ x M x M"^ and A^ > 1, 

< Cmilx - x'|2™ + \y- y'|2"^ + \z - zf"^ 

+ Eiixl" - x,^|2™ + \Yf - y^^i^H). 

Recall that the coefficients $ and / are Lipschitz continuous and the assump- 
tion (E). From Corollary 4.1 and a standard estimate for forward SDEs we 
then get (ii). 

Then, due to Kolmogorov's weak compactness criterion for multi-parameter 
processes the estimates (i) and (ii) of the present lemma and those of Lemma 
4. 1 imply the tightness of the sequence of laws of ^ , A^ > 1 , on C D 

For proving Proposition 4.3 we have still to show the following. 

Lemma 4.5. The finite- dimensional laws of the stochastic fields ^^ = 
(C"^'^,'^^'^,C^'^,C^'^), N > 1, converge weakly to the corresponding finite- 
dimensional laws of the continuous zero-mean Gaussian multi-parameter 
process ^ = {^^^\^^'^\£}-^\£}'^') introduced in Theorem 4-2. 

Proof. Following the idea of the proof of Lemma 4.2 we decompose 
^N ^ (^i,7V^^2,Af^ ^3,iV^^4,iV) j^^^ ^j^g ^^^ ^^ ^^le random fields C'^'^ = (C^'^'^ 
^2,N,i^ ^3,JV,i^^4,7V,i) ande'^'2 = (^^i,N,2^^2,N,2^^3,N,2^^A,N,2^^ where the fields 
^^' '^ and ^^'^'^ have been introduced in the proof of Lemma 4.2, the fields 
,^^'^'^ and ^^^'^'2 are defined in the same way, with b instead of a, and ditto 
for ^3>^'i and ^■^'^'^ where we have $ instead of a and t = T. Finally, in the 
same spirit we define 

1 ^ 
ef'^'^(A) = ^(e'^(/(A, A,)) - E[f{X,At)]), 

^^ k=i 

(t,A)G [0,T] xM'^xMxM'^; 

1 ^ 
ef'^'2(x) = ^=^(e'=(/(A,Af ) -/(A,A,)) -i?[/(A,Af ) -/(A,A,)]) 

for (t. A) G [0, T] X M'^ X M X M'^. 

Since the bounded multi-parameter fields @''{b{-,X),a{-,X),^{-,XT), /(•, 
A)), k>0 are i.i.d., it follows directly from the central limit theorem that 
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the finite-dimensional laws of the random field ^'''^ = {^^' '^, ^2,^,1 ^^3, TV,! ^ 
^4:,N,i^ converge weakly to the corresponding finite-dimensional laws of the 
zero-mean Gaussian multi-parameter process ^ = {C ,C jC j C ) whose 
covariance functions coincides with those of the field (b{-,X), a{-,X), ^(-jXt), 
/(•,A)) [recah the definition of ^ = (^(^), ^(2) ^^(3)^^(4) ^) given in Theorem 4.2]. 

On the other hand, by repeating the argument developed in the proof of 
Lemma 4.2 and recalling the assumption (E) we get from Proposition 3.2 
and Theorem 3.2 that 

^[le^^•^(-)l^] +i^[|e.^'^'^(x)P] +i?[|C^'^'^(x)P] +i?[|C^^'^(x,y,.)|^ 

4 = 1 *J = 1 

+ Var($(x, X^) - $(x, Xt)) 

+ Yav{f{{x,y,z),A^)-f{{x,y,z),At)) 

< C{E[\X^ - Xt\'] + E[\X^ - Xt\^] + E[\Yf - Ytf]) 

< — foT all N> I, {t,x,y,z)e [0,T] x R'^ x M x M*^. 

Consequently, it follows that also the finite-dimensional laws of the random 
fields C^ = {^^'^,f'^, C^'^,^^'^), N>1, converge weakly to the corre- 
sponding finite-dimensional laws of the zero-mean Gaussian field ^ = {S,^^' , 
^{2)^^(3)^^(4)^)^ The proof is complete. D 

As we have already pointed out above. Proposition 4.3 follows directly 
from the Lemmas 4.4 and 4.5. On the other hand. Proposition 4.3 allows to 
adapt the approach for the study of the limit behavior of \/N{X — X) to 
that of the triplet y/N{X^ -X,Y^ - Y, Z^ - Z): using the same notation 
as in the preceding section we remark that, since the random fields ^ , 
N >1, are independent of the driving Brownian motion W , Skorohod's 
representation theorem yields that on an appropriate complete probability 
space {Q',T',P') there exist copies ^'^ , N >l, and ^' of the stochastic 
fields ^ , N > 1, and ^, as well as a d-dimensional Brownian motion W' , 
such that: 

(i) P^,=P^,PI:,^=P^N,N>1; 

(ii) ^'N = (^'i,^^^'2,7V^^/3,^^^/4,W) ^e = (^'^e'^e'^e'^), uniformly on 
the compacts of [0,r] x M"* x M x M*^, P'-a.s.; 

(iii) W is independent of £,' and £,'^ , for all A^ > 1. 

The probability space {0,' ,J-' ,P') is endowed with the filtration F' = {J^l = 
•^r V.^o)te[o,T], where T[, = a{e(x),Cf (x), {s,x) G [0,r] x M, A^ > 1} VAA. 
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Remark that, with respect to this filtration, the process W' is a Brownian 
motion and has the martingale representation property. 

Given ^ , ^' and the Brownian motion W' we redefine the triplets {X ,Y , 
Z^) and {X,Y,Z) on the new probability space: the processes X^ and X 
are redefined as in the preceding section (their redefinitions are denoted 
again by X'^ and X' , resp.), and in the same spirit we introduce the cou- 
ples {Y'^,Z'^), {Y',Z')eSl,{[0,T];R) x L|,([0,r];M'^) as unique solution 
of the backward equations 

yr = [^e'''ix^'')+E'Mx,xni.=x'-) 

+ / (^e"^(Af ) + E'ifiX, Af )]|,=A,.) ds - j^ Zf dW;^, 

tG[0,T] 



and 



Yl = E'mx,X^)],^x^+J^ E'[f{X,A'^)U=^,ds-J^ Z[dWl, 

tG[0,T], 

respectively, where A'^ = (X'^, F'^, Z'^) and A' = {X',Y\Z'). Recah that 
the coefficients ^'^'^(•) and ^s'^(-) are Lipschitz, uniformly with respect to 
(a;,s) G J7 X [0,r], and so are C'^'^(-) and C^^(-) with respect to {J,s) G 
9.' X [0,T]. Thus, since P' o (^'^,VF')"^ = P o {^^ ,W)-^ , N>1, and P' o 
(?', W)-^ = Po{C, W)-^, we can conclude that also P' o (^'^, W' , A'^)-^ = 
Po (^^, W, A^)-\ for all A^ > 1, and P'o (^', W' , A)"! = Po (^, W, A)-^ This 
allows to reduce the study of the limit behavior of \/N{A — A),N > 1, to 
that of the sequence \/]V(A'^ - A'), A^ > 1. 

Proposition 4.4. Let X G 5|,([0,T];M'^) be the unique solution of SDE 
(4.5) and (y,i) G5|,([0,r];M) x L|,([0,r];]R'^) that of the backward equa- 
tion 

Yt = leiX^) + E'[V,^{x,Xlr)l=x^XT 
+ E'[V,>^x,X't)Xt]^=x!,} 

(4.12) +j\^',\K)+E'[Vxf{X,K)]x=A'As 

+ E[Vx'f{X,K)As]x=AOds 
-f ZsdWs, te[0,T],whereA={X,Y,Z). 
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E' 



sup (|^(xf - X',) -Xs\' + \Vn{y:^ - yi) - y. 



r/N 



.se[o,r] 



rT 1 

+ \^{Zf -Z',)-Zs\^ds — ^0 
Jo 

as N ^ +00. 



Proof. We begin by remarking that an argument similar to that of Step 
1 of the proof of Proposition 4.2, combined with estimate (i) of Lemma 4.4, 
allows to show that i'^{X'^) G LP(0', J^^,P') and ^'^(A') G L^,([0,r];M), for 
all p>l. Consequently, since the coefficients V^:^, V^'^, Va/ and V \> f are 
bounded, the above backward SDE admits a unique solution {Y,Z). More- 
over, this solution belongs even to cS|;,([0,T];M) x L^,(fi';L2([0,r];M'^)), for 
anp> 1. 

Prom Step 2 of the proof of Proposition 4.2 we know already that 
(ef'^(Xr))ie[o,T] - (^f (^0)*e[o,T] in LH[0,T] x n',dtdn, as N^+oo. 
The same argument yields that also {^''-^'^ (X!^)) -^ (^'(3)(X^)) in L'^{n',T', 
P'), as N^ +00. Moreover, s £ [0,T], 

^'[ie^(Af)n 

•id'^(A)r 



< E 



sup 

xRxI 



2m- 



1/2 / ^ 1/2 

i?[(i+iAfm 



1 + |A| 
Thus, using Lemma 4.4(iii) and the fact that, for all A^, tti > 1, 

EWA^l^n < CmEWX^l^"' + ly^P'" + iZf |2™1 < Cm, ds-a.e 



Recall that the coefficients of (4.1) and of BSDE (4.7) are bounded, uni- 
formly with respect to A^ > 1; the same holds true for Z (cf. Lemma 4.3), 
we can conclude that, for some constant Cm, 



■fA,N r A /TV M mi 



E'[\C^i^. 



On the other hand, recalling that ^'^ 
C'^)C'^)?'^) uniformly on the compacts of [0, T] x 
A^> 1, 



<Cm, ds-a.e. on [0,T],N>1. 



, P'-a.s., and 



E' 



sup (|xf -xi\' + |yr - Y^\') + r\z','' - z[\^ 

efo.Tl JO 



5e[o,T] 
= E 



sup (|Xf 
.se[o,T] 



\" |2 I IV^ 



Ys?) + 



\Z, 



N 



dt 



Ztl^dt 



< 



C_ 

N 



MFBSDES: A LIMIT APPROACH 39 

(cf. Proposition 3.2 and Theorem 3.2), we see that 

Y: sup \^'^'^iX'f)-C{X:)\ + \e'''{X^'')-e{X!r)\^0 
i=ise[o,T] 

in probabihty P' , and ^^^'^(A',^) — >C^{K), in measure dsdP'. 

Combining this with the uniform square integrabihty of {^'J^ {X'J^),s G 
[0,T],N> 1}, proved above, yields the convergence of ^'*'^(X'^) to C(^') 
in L2([0, T] X n', dtdP'), for i = l,2, that of ^'^'^(X'^) to e'^(X') in L'^{n',r,P') 
and that of ^'^'^(A'^) to ^'^(A') in L2([0,T] x Vt',dtdP'). 

Let us now prove the convergence of ^/N{X'^ - X' ,Y'^ - Y' , Z'^ - Z') 
to A = {X,Y,Z). From Proposition 4.2 we know aheady that 



E' 



sup \VN{X'f - Xi) - Xs\^ 
se[o,t] 



■0 asiV^+oo. 



For verifying the convergence of \/N{Y — Y' , Z — Z') we observe that 

Viv(y/^ - y/) -Yt = il'"" + /2'^ + /3'^ - fiVNiz',"" - z^) - zj dw^, 

te[o,T], 

where 

4'^ = {^'''^{X^^) - ^iX^)} + y[^(e;^'^(Af ) - e'^(A'J) ds, 

/2.^ = VN{E'Mx,X!r'')U^^N - E'Mx,Xlr)Uj,^) 

-iE'[V.Mx,X^)l=x^XT + E'[V.,,<^ix,X^)XTl=x'J, 

/i''^ = ^'^{ViV(i?'[/(A,Af)],=^,.-i?'[/(A,A'J],=A,) 

- {E'[Vxfi\,K)]x=A'As + E'[Vx'fiX,K)As]x=AO}ds, 

t£[0,T]. 

From the convergence of C^^^{X'^) to C'^(X') in L'^{n' ,J^' ,P') and that of 
e'^'^(A'^) to C'^(A') in L2([o,r] X n',dtdP') we obtain that 



S' 



sup |/j 

tG[0,T] 



1,7V|2 



•0 asiV^+cx). 



In analogy to Step 3 of the proof of Proposition 4.2 we get that, for some 
real sequence < pj\f\0 (X — > +oo) , 

2i 



E'[\l''^f]<PN + CE' 



Vn{X^ -X't)-Xt 
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<PN + CE' 



sup \Vn{xI^^ -X'^)-Xt\ 

.te[o,T] 



as A^ ^ +00. 



It remains to estimate I^' . For this, thanks to the continuous difFerentia- 
bihty of /, with the notation A'/''^ := A'^ + 7 (A'/ - A'J we get 



E' 



sup |/f'^|2 

.se[t,T] 



<CE' 



T 



E'[Vxf{Xi,AT'') 

- Vxf{X2,K)]\x,=AT'-'M=A', d^ X ^s 



]lX,=A'f'\X,=A'J^ 



ds 



+ /'i5;'[{Vv/(Ai,Af'^) 

JO 

-Vv/(A2,A'J}A, 
+ C{T-t)E' 
X j^ \E'[Vxf{\M^^^)]x=A'S^{M^ - K) - A.) 

+ E'[Vx'fi\Af'')i^iK''-K)-h)]x=A'fds , 

t€[0,T], 

and taking into account that the first-order derivatives of / are bounded 
and continuous, it fohows from the convergence of A •'^ — > A' (A^ -^ +cxd) in 
i3|/, for all 7 S [0, 1], that for some sequence < pat \ 0, 



E' 



sup |/;, 



3,N\2 



<PN + C{T - t) 1^ E'[\VN{Af - A'J - As\^]ds, 



te[o,r],iV>i. 



Thus, for some sufficiently small 5 > which depends only on the Lipschitz 
constant of /, and for some sequence < p^ \ 0, we have for all t £ [T — 

E'[\^{Y:''-Y:)-Ytf] + ^E' 



N{Z'f - Z',) - Zs\ ds 



<PN + CE' 



f \^{X'f -X'^)-Xsfds 



+ CE' 
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YD-YJ^ds 



'N{Y^^ 

and from Proposition 4.2 and Gronwall's inequality we get 

E'[\^/N{Y^^ - y/) - fi I ^] ^ for all te[T-6,T] 
and 



E' 



T~S 



N{z: 



/N 



ZA — ZJ ds 







as A^ ^ +00. 



Having got the convergence of y/N{Y^^g - Y^_g) - Yts in L'^{n' ,T' ,P') 
we can repeat the above convergence argument on [T — 25, T — 5] by replacing 
P'^ by VNiY^'Ls-YT^s) - Yt^s. This yields E'[\VN{Y/^ - Y/) - Yt\^] ^ 0, 
for all t £ [T - 26,T - 6], and E'[J^Z2S l^l^^ - Z'^) - Zs\^ ds] ^ 0. Iter- 
ating this argument we get the convergence over the whole interval [0,r]. 
Finally, from the equation 

^/]v(y/^ - y/) -Y, = //'^ + /2'^ + 1^^ - l^iVNizf - zi) - zj dwi, 

t€ [0,r], we get 
E' 



sup |Viv(y/^-y/)- 

e[o,T] 


-Ytf 












<4LB' 
V 


1 rl.^|2 
sup J( 

.tG[0,T] 


+ ^'[|J2,7V|2] 








+ E' 


' sup |/^^|2' 
-te[o,T] 


+ i5;' 


' r\^{z[^-z[)- 

Jo 


-Zt\^dt 


) 


<PN + 


C7. 


E' 


Jo 


(A 


/N 
s 


-A 


's)-h?ds 


— >0 


as N 


-^ +00. 



The proof is complete. D 

Finally, it remains to give the proof of Theorem 4.2. 

Proof of Theorem 4.2. In analogy to the existence and uniqueness 
of the triplet A = (X, Y, Z) £ fi|, we can prove that of A = (X, F, Z) G B^. 
Since W is independent of ^ as well as C^ ,N > 1, and, on the other hand, 
also W' is independent of ^' and ^'^ ,N > 1, it follows that: 



(i) P o {W,C^)~^ = P' o {W' ,C ) , for alliV>l; 
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(ii) po{w,i)-^=p'o{w\er\ 

We notice that A'^ is the unique solution of an forward-backward SDE 
governed by the W' and ^'^, and A^ is the unique solution of the same 
system of equations but driven by W and ^ instead of W' and ^ . On the 
other hand, A' is the unique solution of the system (4.8)-(4.9) with W' at 
the place of W . Thus, 

(iii) Po (l^,e^,A,A^)-i = P'o (T^',^'^,A',A'^)-i, for all iV> 1. 

Concerning the process A, it is the unique solution of the same equation 
as A, but with ^ and W replaced by ^' and W . This, together with (ii), 
yields: 

(iv) Po(l^,e,A,A)-i = P'o(M/',e',A',A)-i. 

From (hi) and (iv) we see, in particular, that -P^/^vca'V-A) = ^'./ni/^in _i^,y 
N >1, and P^ = P^. For completing the proof it suffices now to remark that 

the convergence of y/N{A'^ - A') to A in cS|, ( [0, T] ; R"') implies that of their 
laws on C([0,r];]R'^) x C([0,r];M) x L^{[0,T];R'^). The proof is complete. 
D 
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